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Both the basic electronic structure of tetraborides, and the changes across the lanthanide series 
in _RB4 [R = rare earth) compounds, are studied using the correlated band theory LDA+U method 
in the all-electron Full Potential Local Orbital (FPLO) code. A set of boron bonding bands can 
be identified that are well separated from the antibonding bands. Separately, the "dimer B" 2pz 
orbital is non-bonding [mz. graphite and MgB2), and mixes strongly with the metal 4d or 5ci states 
that form the conduction states. The bonding bands are not entirely filled even for the trivalent 
compounds (thus the cation d bands have some filling), which accounts for the lack of stability of 
this structure when the cations are divalent (with more bonding states unfilled). The trends in the 
mean 4/ level for both majority and minority, and occupied and unoccupied, states are presented 
and interpreted. 



I. BACKGROUND AND MOTIVATION 

The tendency of the metalloid boron to form clus- 
ters has led to widespread study of the properties 
of condensed boron. Of the many classes of com- 
pounds that B forms, B-rich metal borides include 
classes with very important, and intensely studied, 
properties. One example is MgB2, which is the pre- 
mier phonon-coupled superconductor^ (at 40 K). Al- 
though this structural class includes several transi- 
tion metal borides and other simple metal borides 
(such as LaB2), MgB2 is unique in this single- 
member class of quasi-two-dimensional s-p metal 
with very high superconducting transition temper- 
ature due to strong covalent B-B bonds that are 
driven metallic^ by the crystal structure and chem- 
istry. 

Another class that has received great atten- 
tion is the hexaborides MBg formed from vertex- 
linked Bg octahedra that enclose the metal ion 
in the cubic interstitial site. This class includes 
the divalent metals (M=Ca, Sr, Ba) that are 
small gap semiconductors ^liiS ,6 ,7,8 ,9 , lo, 1 1,12 r^j^^ 

bility of this structure was understood decades ago, 
when cluster studies established^!^ that the bond- 
ing states of linked Bg clusters are filled by 20 
electrons, which requires two per Bg unit in addi- 
tion to the B valence electrons. There are many 
trivalent hexaborides as well, including lanthanide 
members which have very peculiar properties: un- 
usual magnetic ordering, heavy fermion formation, 
and s up er conduct i vityi^i^ i 11,13,14,15,16 ^p-^^o monova- 
lent members, NaBg [l3| and KBg have been 
reported. 

Yet another class that has been known for decades 



is the metal (mostly rare earths) tetraboride -RB4 
family, which is richer both structurally and elec- 
tronically and for which considerable data is avail- 
able (see: for several BB^, Rcfs. [19,20,21,22]; YB4, 
Refs. l^llllillilll; LaB4, Ref. [28|]; CeB4, Rcfs. 

fllsOml; NdB4, Ref. ^\ GdB4, Refs. i33,M 
36 37 SSh TbB4^efs. l3l|40l4l||42|4iraiDyB4, 
Refs. |45ii4ij43|4i|45,'5S5; ErB4, Refs. |4l,l5T,'5^ ) . 
Yttrium and all the lanthanides except Eu and Pm 
form isostructural metallic tetraborides i?B4 with 
space group P4/mbm (#127), described below and 
pictured in Fig. [TJ Presumably Eu is not stable in 
the tetraboride structure because of its preference 
for the divalent configuration in such compounds. 
The Sr and Ba tetraborides also are not reported. A 
"calcium tetraboride" with formula Ca(Bi_a;Ca;)4, 
X « 0.05 was reported^ recently. 

These rare-earth tetraborides exhibit an unusual 
assortment of magnetic properties. While CeB4 
and YbB4 (/^ and f"^^ respectively) don't order 
and PrB4 orders ferromagnetically at Tc=25 K,'^- 
all of the others (i?=Nd, Sm, Gd, Tb, Dy, Ho, 
Er, Tm) order antiferromagnetically, with Neel tem- 
perature T^v (see Table I) spanning the range 7- 
44 K. A noteworthy peculiarity is that T^r doesn't 
obey de Gennes' scaling law, which says that the 
magnetic transition temperature is proportional to 
{gj — 1)^J(J -|- 1) across an isostructural series 
where the rare-earth atom is the only magnetic 
component4^i^ (here J is the Hund's rule total 
angular momentum index, gj is the corresponding 
Lande ij-factor.) In the rare earth nickel borocar- 
bide series, for example, de Gennes scaling is obeyed 
faithfuUy.!^^ This lack of scaling indicates that mag- 
netic coupling varies across the series, rather than 
following a simple RKKY-like behavior with a fixed 




FIG. 1: (color online) Structure of KB4, viewed from 
along the c direction. The large metal ion spheres (red) 
lie in z=0 plane. Apical Bl atoms (small black) lie in z 
~ 0.2 and z ~ 0.8 planes. Lightly shaded (yellow) dimer 
B2 and equatorial B3 (dark, blue) atoms lie in z—0.5 
plane. The sublattice of R ions is such that each one is 
a member of two differently oriented R4 squares, and of 
three R3 triangles. 



Fermi surface. 

Both the ferromagnetic member PrB4 and an- 
tiferromagnetic ones RB4 show strong magnetic 
anisotropy. For ferromagnetic PrB4 the c axis is the 
easy axis. The situation is more complicated for the 
antiferromagnetic compounds, which display vary- 
ing orientations of their moments below Tn, and 
some have multiple phase transitions. GdB4 and 
ErB4 have only one second order phase transition, 
while both TbB4 and DyB4 have consecutive second 
order phase transitions at distinct temperatures. A 
yet different behavior is shown by H0B4 and TmB4, 
which have a second order phase transition followed 
by a first order phase transition at lower tempera- 
ture. The magnetic ordering temperatures, primary 
spin orientations, and experimental and theoretical 
effective (Curie- Weiss) magnetic moments have been 
collected in Table L 

The variety of behavior displayed by these tetra- 
borides suggests a sensitivity to details of the un- 
derlying electronic structure. Unlike the intense 
scrutiny that the tetraborides have attracted, there 
has been no thorough study of the tetraboride elec- 
tronic structure, which contains a new structural el- 
ement (the "boron dimer" ) and an apical boron that 
is inequivalent to the equatorial boron in the octa- 
hedron. We provide here a detailed analysis, and in 
addition we provide an initial look into the trends to 
be expected in the 4/ shells of the rare earth ions. 



TABLE I: Data on magnetic ordering in the _RB4 
compounds.— i^ii^iS The columns provide the experi- 
mental ordering temperature(s) Tmag, the ordering tem- 
perature Tth predicted by de Gennes law (relative to the 
forced agreement for the GdB4 compound), the orienta- 
tion of the moments, and the measured ordered moment 
compared to the theoretical Hund's rule moment (/xs). 







Tth (K) 


direction 


M(exp) 


M(th) 


PrB4 


24 


2.1 


II c 


3.20 


3.58 


SmB4 


26 


12 






0.84 


GdB4 


42 


42 


_L c 


7.81 


7.94 


TbB4 


44, 24 


28 


_L c 


9.64 


9.72 


DyB4 


20.3, 12.7 


19 


II c 


10.44 


10.63 


H0B4 


7.1, 5.7(lst) 


12 


II c 


10.4 


10.6 


ErB4 


15.4 


7 


II c 


9.29 


9.60 


TmB4 


11.7, 9.7(lst) 


3 


_L c 


7.35 


7.56 



II. CRYSTAL STRUCTURE 

The full i?B4 structure was first reported by Za- 
Ikin and Templeton^ for the Ce, Th, and U mem- 
bers. These tetraborides crystallize at room tem- 
perature in the tetragonal space group P4-/mbm, 
D\f^ with four formula units occupying the positions 
listed in Table HT] The lattice constants for the re- 
ported rare earth tetraborides are presented in Table 

[ml 

The Bl and B3 atoms form Be octahedra (apical 
and equatorial vertices, respectively) that are con- 
nected by B2 dimers in the z=l/2 plane. The Bg 
octahedra, which are arrayed in centered fashion in 
the x-y plane within the cell, are flattened some- 
what, with distances from the center being 1.20 A 
along the c axis and 1.29 A in the x-y plane (tak- 
ing GdB4 as an example). Each B2 atom is bonded 
to two Bl atoms in separate octahedra and to one 
other B2 atom. A suggestive form for the chemical 
formula then is [i?2B2B6]2- The rare-earth atoms 
lie in the large interstitial holes in the 2=0 plane, 
and form a 2D array that can be regarded as fused 
squares and rhombuses. '^^ 

The R site symmetry is mm. The symmetry of 
an R site is important for the properties of the com- 
pounds, as it dictates the crystal field splitting of the 
ion with total angular momentum J = L + S and 
thereby the resulting magnetic state at low temper- 
ature. The R ion is coordinated by seven B atoms in 
planes both above and below, three of them being 
dimer B2 atoms (two 2.88 A distant and one at a 
distance of 3.08 A) and four of them equatorial B3 
atoms (two each at distances of 2.76 A and 2.84 A). 
Within the unit cell the four R sites form a square 
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FIG. 2: Plot of experimental lattice constants of RB4, vs 
position in the Periodic Table (atomic number), showing 
a lanthanide contraction of about 5% for a, 3% for c. The 
smooth lines show a quadratic fit to the data. 



of side d = 0.518a = 3.70A, oriented at about 15° 
with respect to the square sublattice of Be octahe- 
dra. The (low) site symmetries of the apical Bl, 
dimer B2, and equatorial B3 atoms are 4, mm, m, 
respectively. 

The reported lattice constants for the lanthanides 
are plotted in Fig. O It is evident that most fall 
on smooth lines reflecting the lanthanide contrac- 
tion in this system. The behavior is representative 
of trivalent behavior, from La through to Lu. The 
big exception is Ce, which has smaller volume sug- 
gesting that, rather than being simple trivalent, the 
4/ electron is participating in bonding. Pm with all 
unstable isotopes has not been reported. EUB4 also 
has not been reported; Eu typically prefers the diva- 
lent state (due to the gain in energy of the half-filled 
4/ shell) so it is not surprising that it is different. 
However, some divalent tetraborides do form in this 
structure {e.g. CaB4, see Sec. IV) so it cannot be 
concluded that EUB4 is unstable simply on the ba- 
sis of divalency. Finally, the small deviation of Yb 
from the smooth curves suggest it maybe be mixed 
or intermediate valent (although close to trivalent). 
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TABLE II: Site designations, symmetries, and atomic 
positions of the atoms in the _RB4 crystal. 

R 4g mm (x, i+x, 0) 
Bl 4e 4 (0, 0, z) 
B2 4h mm (x, i+x, i) 
B3 8j m (x, y, i) 



TABLE III: Tabulation of the lattice constants and inter- 
nal structural parameters used in our calculations. Con- 
sidering the extreme regularity of the internal coordi- 
nates through this system, the irregularity in zbi for Dy 
should be treated with skepticism. 



R 


a(A) 


c(^) 


Xr 


Zbi 


Xb2 


XB-i 




Ref. 


Y 


7.111 


4.017 


0.318 


0.203 


0.087 


0.176 


0.039 


23 


La 


7.324 


4.181 


0.317 


0.209 


0.088 


0.174 


0.039 


22.28 


Ce 


7.208 


4.091 


0.318 


0.203 


0.087 


0.176 


0.039 


22.30 


Pr 


7.235 


4.116 


0.318 


0.203 


0.087 


0.176 


0.039 


21 


Nd 


7.220 


4.102 


0.318 


0.203 


0.087 


0.176 


0.039 


22^ 


Pm 


7.193 


4.082 


0.318 


0.203 


0.087 


0.176 


0.039 




Sm 


7.179 


4.067 


0.318 


0.203 


0.087 


0.176 


0.039 


21 


Eu 


7.162 


4.057 


0.318 


0.203 


0.087 


0.176 


0.039 




Gd 


7.146 


4.048 


0.317 


0.203 


0.087 


0.176 


0.038 


35 


Tb 


7.120 


4.042 


0.317 


0.202 


0.087 


0.176 


0.039 


41.43 


Dy 


7.097 


4.016 


0.319 


0.196 


0.086 


0.175 


0.039 


21,51 


Ho 


7.085 


4.004 


0.318 


0.203 


0.087 


0.176 


0.039 


21 


Er 


7.071 


4.000 


0.318 


0.203 


0.086 


0.177 


0.038 


43.51 


Tm 


7.057 


3.987 


0.318 


0.203 


0.087 


0.176 


0.039 


22 


Yb 


7.064 


3.989 


0.318 


0.203 


0.087 


0.176 


0.039 


22 


Lu 


7.036 


3.974 


0.318 


0.203 


0.087 


0.176 


0.039 


22 



III. CALCULATIONAL METHODS 

The fuU potential local orbital (FPLO) code^^ 
(version 5.18) was used in our calculations. 
Both LDA (PW92 of Perdew and Wang") and 
LDA-I-U (using the atomic limit functional) 
are used. We used a k mesh of 12^ in the 
full Brillouin zone. For the density of states 
(DOS) plot and Fermi surface plot, we used a 
k mesh of 24^ for more precision. The basis 
set was Is2s2p3s3p3d4s4p::(4d4f5s5p)/6s6p5d-h 
for all metal elements (except 

Y(ls2s2p3s3p3d::(4s4p)/5s5p4d-|-) and 
Ca(ls2s2p::(3s3p)/4s4p3d-|-)). For boron atoms we 
used the basis ::ls/(2s2p3d)-|-. 

In the LDA-f U calculations we used values typ- 
ical for 4/ atoms U = 8 eV and J = 1 eV (cor- 
responding to Slater integrals Fi=8.00, F2=11.83, 
F4=8.14, F6=5.86) throughout all calculations. The 
high symmetry points in the tetragonal zone are 

r=(o,o,o), X = (f,o,o), M = (f,f,o), z = 
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IV. GENERAL ELECTRONIC STRUCTURE 



YB^ band structure 
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The valence-conduction band structure of YB4 
(where there are no 4/ bands) is shown in Fig. [3] 
For LaB4, which differs in volume and conduction d 
level position, the bands are very similar, with only 
slightly differing Fermi level crossings along the M-F 
direction. The occupied valence bandwidth is 11 eV 
(not all bands are shown in this figure). One strik- 
ing feature of the bands is the broad gap of more 
than 3 eV along the top (and bottom) edges R-A-R 
of the Brillouin zone. Bands along these lines stick 
together in pairs due to the non-symmorphic space 
group, and nearly all bands disperse very weakly 
with (or ky) along these edges. This gap closes 
along the k^ = tt/c plane of the zone only for small 
in-plane components of the wavevectors. It is such 
gaps enclosing that often account for the stabil- 
ity of a crystal structure, and the stability of boride 
structures, including this one, has been a topic of 
interest for decadesi^i^ i^^'^" 

The band structure of a divalent cation member 
(CaB4) is also included in Fig. [3] for comparison. 
The largest difference is the band filling, as expected, 
although some band positions differ in important 
ways near the Fermi level. Still the 3d bands of 
Ca are not quite empty, as a band with substantial 
3d character lies at Ej? at R and is below F,p all 
along the R-A line. CaB4 can be fairly character- 
ized, though, as having nearly filled bonding B 2p 
bands and nearly empty Ca 3d bands. 



A. Bonding and Antibonding Bands 



FIG. 3: Band structure of YB4 (top panel) and CaB4 
(lower panel) within 6 eV of the Fermi level along high 
symmetry directions, showing the gap that opens up 
around E_f (taken as the zero of energy) throughout 
much of the top and bottom portions of the tetrago- 
nal Brillouin zone. Notice the lack of dispersion along 
the upper and lower zone edges R-A-R (kz—ir/c, and 
either or ky is n/a). Note also that, due to the non- 
symmorphic space group, bands stick together in pairs 
along X-M (the zone 'equator') and along R-A (top and 
bottom zone edges). 



(0,0,f),i?=(^,0,f),andA=(f,f,^). 



As mentioned in the Introduction, the stability of 
the hexaborides is understood in terms of ten bond- 
ing molecular orbitals of the Bg octahedron. This 
octahedron occurs also in these tetraborides, along 
with one additional B2 dimer that is bonded only 
in the layer (sp^). Lipscomb and Britton^i^ started 
from this point, and argued that each of the B2 
atoms in a dimer forms single bonds with two B3 
atoms but a double bond with its dimer neighbor, 
so each B2 atom needs four electrons. The total 
of 20-1-8 electrons for each set of 6+2 boron atoms 
leaves a deficit of four electrons, or a deficit of 8 elec- 
trons in the cell. This amount of charge can be sup- 
plied by four divalent cations, with CaB4 as an ex- 
ample. Most tetraborides contain trivalent cations, 
however, so this is an issue worth analyzing. 

An empirical extended Hiickel band structure 
study^ for CaB4 indeed gave a gap, albeit a very 
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FIG. 4: Projected density of states per atom of each of 
the B atoms for YB4. The curves are shifted to enable 
easier identification of the differences. The B 2p bonding- 
antibonding gap can be identified as roughly from -1 eV 
to 4-5 eV. 




FIG. 5: Enlargement of the partial densities of states of 
Y 4d and B 2p states (per atom) near the Fermi level. 
The states at the Fermi level, and even for almost 2 eV 
below, have strong 4d character. The apical B2 character 
is considerably larger than that of Bl or B3 in the two 
peaks below Ef, but is only marginally larger exactly at 
Ef . 



narrow one. The Hiickel method can be very instruc- 
tive but is not as accurate as self-consistent den- 
sity functional methods. Our FPLO calculation on 
CaB4, shown in Fig. 3, gives a metallic band struc- 
ture. However, the 'valence' (occupied) and 'conduc- 
tion' (unoccupied) bands in the bands (Hiickel, and 
also FPLO) are readily identified, and it clear that 
there are disjoint sets of bands with different charac- 
ters. There are the boron bonding bands (at Fip and 
below) that can be clearly distinguished from con- 
duction bands at and above Fjp- These conduction 
bands are primarily metal d bands (with an inter- 
spersed nonbonding B2 pz band, see below). If they 
were ^0.5 cV higher it would result in an insulat- 
ing band structure in CaB4. The boron antibond- 
ing bands lie higher, above 5 eV at least and mix 
strongly with the metal d bands. 

The separation into bonding and antibonding B 
2p bands agrees (almost) with the ideas of Lipscomb 
and Britton, and confirms their counting arguments. 
However, the existence of numerous R^^^4 com- 
pounds and only one divalent member shows that 
the extra electron is not a destabilizing influence, 
while it increases the conduction electron density 
(hence, the conductivity, and magnetic coupling). 

In covalently bonded materials it is common to be 
able to identify the distinction between the bonding 
bands and the antibonding bands. In covalent semi- 
conductors, for example, they lie respectively below 
and above the band gap, an absolutely clean separa- 
tion. In the KBi system the d bands lie within the 
corresponding bonding-antibonding gap and compli- 
cate the picture. Analysis of the orbital-projected 
bands clarify this aspect. The Bl and B3 atoms, 
being engaged in three-dimensional bonding (within 
an octahedron and to another unit [octahedron or 
dimer]), have a clear bonding-antibonding splitting 
of a few eV (beginning just below E^). Likewise, the 
dimer B2 Px,Py states display a similar splitting. 

The B2 pz orbital is quite different. As is the 
case in MgB2 (whose planar structure is similar to 
the local arrangement of a B2 atom), pz bands ex- 
tend continuously through the gap in the B bond- 
ing/antibonding bands, and mix fairly strongly with 
the rare earth d states in that region. There is con- 
siderable B2 pz character in the bands near (both 
below and above) Ei^- at the zone edge M point, as 
well as the Y 4c? character that is evident in Fig. [3l 
So while there is some Bl and B3 character in the 
rare earth metal d bands that lie within the boron 
bonding-antibonding gap, the amount of B2 pz char- 
acter is the primary type of B participation in these 
bands that provide conduction and magnetic cou- 



pling. 



B. Pseudogap in the Density of States 

From the projected DOS of the three types of B 
atoms of YB4 (see Fig. |4]), one can detect only rela- 
tively small differences in the distribution of Bl, B2, 
and B3 character arising from their differing envi- 
ronments. First, note that in the DOS of Bl and B3 
there is a peak around -15 eV, while there is no such 
peak for B2. This peak arises from the overlap of 2s 
and 2pa- states of each of the boron atoms forming 
the Bg octahedra (Bl and B3); the 2s character is 
about three times as large as the 2p^ character, and 
the remaining 2s character is mixed into the lower 
2p bands. This state is a well localized Bg cluster 
orbital, and there are two such orbitals (octahedral 
clusters) per cell. The bridging B2 atoms do not 
participate in any such bound state. 

Another difference in characters of the B sites is 
that, in the region below but within 2 eV of the 
Fermi level, the DOS of the dimer B2 atom is signif- 
icantly larger than that of Bl and B3 atoms, as can 
be seen in Fig. [51 Together with plots showing the 
band character (not shown), this difference reflects 
the fact that all of the 2p orbitals of Bl and B3 (octa- 
hedron) atoms are incorporated into bonding (filled) 
and antibonding (empty) bands. The distinct char- 
acteristic of the B2 pz state was discussed in the 
previous subsection. All B 2p states do hybridize to 
some degree with the metal d bands, however, and 
all B atoms have some contribution at the Fermi 
level. 

The full Y 4d DOS (not shown) establishes that 
these bands are centered about 4 eV above Ep, with 
a 'bandwidth' (full width at half maximum) of 6-7 
eV (a 'full bandwidth' would be somewhat larger). 
The largest Y character near Ep along symmetry 
lines is 4d(2:^ ~ y^)i primarily in the bands dispers- 
ing up from -0.5 eV at Z toward F. The flat bands 
around -1 eV along T — X — M — T are strongly 
4d(z^) character, indicative of a nonbonding, almost 
localized state in the x-y plane. Note that these 
bands disperse strongly upward along {0,0, kz) and 
lie 3-4 eV above Ej? in the — tt/c plane. Thus 
the 4(i(z^) orbitals form two nearly separate one- 
dimensional bands along k^, and give rise to flat 
parts of some Fermi surfaces (see following subsec- 
tion). These bands can be modeled by a tight- 
binding band —tddaCOskzC with hopping amplitude 
tdda ~ 1 eV. Most of the Ad{xz),Ad{yz) character 
and 4:d{xy) character lies above Ep, and is centered 





FIG. 6: (color online) Fermi surfaces of YB4. Light 
(yellow) surfaces enclose holes, dark (red) surfaces en- 
close electrons. The wide gap the throughout the top 
and bottom edges of the zone account for the lack of 
Fermi surfaces there. 



3-4 eV above Ep. The B2 2pz state mixes primarily 
with Y 'idxzi '^dyz near the M point (near Ej? and 
above) . The B2 2pz orbitals are shifted up somewhat 
with respect to the 2px, 2py states by the ligand field 
effects (there is a bonding interaction within the x-y 
plane only). 



C. Fermi Surface 

The Fermi surfaces of YB4, shown in Fig. [6l will 
be representative of those of the trivalent i?B4 com- 
pounds although small differences may occur due 
to element-specific chemistry of trivalent rare earths 
and due to the lanthanide contraction. The large gap 
along the R-A-R edges precludes any FS on or near 
most of the fc, = - face. The Fermi surfaces can be 

^ c 

pictured as follows. Square hole pyramids with only 
slightly rounded vertices lie midway along the T — Z 
line, and similar nested electron pyramids lie along 
the M — A line near the M point. A pointed ellipsoid 
oriented along kz sits at the Z point. Surrounding F 
is lens-type electron surface joined to pointed ellip- 
soids along the (110) directions. Finally, there are 
two "tortoise shell" shaped hole surfaces within the 
zone, centered along the T — Z lines. 

These surfaces, and the small variation through 
the lanthanide series, is surely relevant to the vary- 
ing magnetic behavior observed in i?B4 compounds. 
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FIG. 7: The full valence band structure of DyB4, and 
up to 5 eV in the conduction bands. This plot is for 
ferromagnetic alignment of the spin moments, with the 
solid bands being majority and the lighter, dashed lines 
showing the minority bands. The flat bands in the - 
4.5 eV to -11 eV are 4/ eigenvalues as described by the 
LDA+U method. 



There are nesting possibilities between the bases of 
the square pyramids, for example, which will ap- 
pear as RKKY coupling as the associated nesting 
vectors. The ellipsoidal attachments on the zone- 
centered lens surface may provide some weak nest- 
ing. 



V. THE LANTHANIDE SERIES 

Any effective one-electron treatment of the elec- 
tronic structure of 4/ electron systems faces seri- 
ous challenges. The root of the difficulty is that the 
ground state of an open 4/ shell has intrinsic many- 
body character, being characterized by the spin S 
and angular momentum L of all of the 4/ electrons, 
and the resulting total angular momentum J, fol- 
lowing Hund's rules. Although it is possible to delve 
into the extent to which the LDA+U method can re- 
produce the z-components of such configurations,^^ 
that is not the intention here. LDA-I-U reliably gets 
the high spin aspect, which contains much of the 
physics that determines relative 4/ level positions 
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FIG. 8: Band structure of DyB4 on a fine scale around 
the Fermi energy, see Fig. [7] The exchange splitting 
(between solid and dashed bands) gives a direct measure 
of the coupling between the polarized Dy ion and the 
itinerant bands (see text). 

and hence trends across the series. There is recent 
evidence from calculations on rare earth nitrides''— 
that, if spin-orbit coupling is neglected and the sym- 
metry is lowered appropriately, the high orbital mo- 
ment (Hund's second rule) can usually be repro- 
duced. The exceptions are the usual difficult (and 
interesting) cases of Eu and Yb. 

The Hund's rule ground state of the ion often 
breaks the local symmetry of the site, and if one is 
exploring that aspect the site symmetry should be 
allowed to be lower than the crystalline symmetry. 
As stated, we are not interested here in those details. 
In the calculations reported here, the crystal symme- 
try is retained. The site symmetry of the lanthanide 
ion is already low (mm), reflected in its 14- fold coor- 
dination with B atoms. In addition, spin-orbit cou- 
pling has not been included. 



A. Band structure 

Most of the RB4 lanthanide tetraborides follow 
the usual trivalent nature of these ions, and the itin- 
erant parts of their band structures are very similar 
to those of YB4 and LaB4. The exceptions are R = 
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FIG. 9: Calculated mean 4/ eigenvalue position (sym- 
bols connected by lines) with respect to Ef, and the 
spread in eigenvalues, of KB4 compounds. The smooth 
behavior from Pr to Tm (except for Eu) reflects the 
common trivalent state of these ions. Eu and Yb are 
calculated to be divalent and deviate strongly from the 
trivalent trend. Ce has a higher valence than three, ac- 
counting for its deviation from the trivalent trend. 

Eu and Yb, which tend to be divalent to achieve a 
half-filled or filled shell, respectively. 

By way of illustration of the complexity of the full 
RB4 bands, the full band structure of DyB4 is pre- 
sented in Fig. [7] for ferromagnetic ordering. The 
4/ bands themselves can be identified by their flat 
(weakly hybridizing) nature. An enlarged picture 
of the bands within 1 eV of Ep is given in Fig. 8. 
The splitting of the majority and minority itinerant 
bands provide a direct measure of the Kondo cou- 
pling of the 4/ moment to the band states. Note 
that the sign of this splitting can vary from band to 
band. 

Figure 8 suggests that the Fermi surfaces will be 
different in the magnetic tetraborides (compared to 
YB4) in specific ways. For Dy, the F-centered sur- 
face splits almost imperceptibly. The surfaces that 
cross the F-Z line also are relatively unaffected by 
exchange splitting. At the M point, however, a new 
surface appears due to the magnetism: an electron 
surface of minority spin. For this band, the polar- 
ization is opposite to the direction of the Dy spins. 
This figure is specifically for ferromagnetic align- 
ment, while DyB4 actually orders antiferromagneti- 



cally (see Sec. I). 

B. Position of 4/ Levels 

The mean position of 4/ levels is displayed in Fig. 
[HI separated into occupied and unoccupied, and ma- 
jority and minority, and trends are more meaningful 
than absolute energies. Simple ferromagnetic align- 
ment is used here, in order to follow the chemical 
trends in the simplest manner. For the occupied 
majority states, the 4/ level drops rapidly from Pr 
(-3 eV) to Sm (-7 eV), then becomes almost flat for 
Gd-Tm (around -8 eV). For the unoccupied minority 
states, the mean 4/ level drops almost linearly from 
Pr {+5 eV) to Er {+2 eV), and for Tm the 4/ level 
is very close to E^. The unoccupied majority levels, 
which become occupied minority levels beyond the 
middle of the series, drop more steeply, with slope 
almost -1 eV per unit increase in nuclear charge. 

There are the usual exceptions to these overall 
trends. Ce is very different, indicating that it is very 
untypical (the calculational result is tetravalent and 
nonmagnetic). Both Eu and Yb are divalent in the 
calculation; an 'extra' 4/ state is occupied so their 
mean 4/ level position is 6 eV (8 eV for Yb) higher 
than the trivalent line. 

The spread in 4/ eigenvalues is also displayed in 
Fig. [31 This spread is sensitive to the specific con- 
figuration that is obtained, and also has no direct 
relation to spectroscopic data, although it does re- 
flect some of the internal potential shifts occurring 
in the LDA-I-U method. The distinctive features are 
unusually large spread for the occupied majority lev- 
els in Dy (two electrons past half-filled shell), and for 
the unoccupied minority (and also unoccupied ma- 
jority) levels in Pr (two electrons above the empty 
shefl). 

VI. SUMMARY 

In this paper we have provided an analysis of the 
electronic structure of trivalent tetraborides, using 
YB4 as the prototype, and compared this with a 
divalent member CaB4. In agreement with earlier 
observations on the likely bonding orbitals in the B 
atoms, it is found that bonding states are (nearly) 
filled and antibonding states are empty. The states 
at the Fermi level in the trivalent compounds are 
a combination of the (dimer) B2 pz nonbonding 
orbitals whose bands pass through the bonding- 
antibonding gap, and the cation d orbitals. Since 



the extra electron in the trivalent compounds does 
not go into an antibonding state, there is no signifi- 
cant destabihzation of the crystal structure. 

The trends in the energy positions of the 4/ 
states in the rare earth tetraborides has been found 
to be consistent with expectations based on other 
rare earth systems, as is the fact that Eu and Yb 
tend to be divalent rather than trivalent. Investiga- 
tions of the magnetic behavior of rare earth tetra- 
borides will require individual study. Nearest neigh- 
bor magnetic interactions may involve a combina- 
tion of 4/ — 4(i — 2pz — 4(i — 4/ interactions, and 
longer range RKKY interactions that may bring in 
the Fermi surface geometry. Another possible cou- 
pling path is the direct 4/ — 2pz — 4/ path. The 
coupling is likely to be even more complicated than 
in the rocksalt EuO and Eu chalcogenides, where 
competition between direct and indirect magnetic 
coupling paths has received recent attentioufS^ The 
tetraboride structure is fascinating in several re- 
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spects. A relevant one, if coupling does proceed 
directly via 4/ — 2pz — 4:f, is that the (dimer) B2 
atom coordinates with three neighboring rare earths 
ions, which will introduce frustration when the in- 
teraction has antiferromagnetic sign. 
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